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ABSTRACT 

We discuss our studies of polarized parton distributions which are related 
to the RHIC-Spin project. First, the parametrization of unpolarized parton 
distributions is explained as an introduction to general audience. Second, 
activities of the RHIC-Spin-J working group are reported. In particular, pre- 
liminary results are shown on the parametrization of polarized parton distri- 
butions. Third, we discuss numerical solution of the DGLAP Q 2 evolution 
equations for longitudinally polarized and transversity distributions. An effi- 
cient and accurate computer program is important for analyzing future RHIC 
experimental data. 

1 Introduction 

Internal spin structure of the proton has been investigated extensively for the last ten years. 
Although the proton spin is no more considered at the stage of "crisis", the clear explanation 
is still far from the final one. In particular, there is essentially no information on sea-quark 
and gluon polarizations. Furthermore, parton's angular momenta may have influence on the 
issue. In order to clarify these problems, the RHIC-Spin project was proposed |]. Obtained 
data at the RHIC facility should be able to pin down how the proton spin consists of internal 
constituent spins. 

The RHIC facility will provide us fruitful information on the proton spin structure. 
In order to understand the situation of proton-spin studies and also to explore possible 
RHIC-Spin physics, we established a Japanese working group on parametrization of polarized 
parton distributions 0. For the time being, the purpose is to obtain the optimum parton 
distributions for explaining existing polarized experimental data. Our study is still at the 
preliminary stage M. The parametrization of unpolarized parton distributions has been 
well studied since 1980's. It is also fortunate that abundant experimental data are available 
through various reactions in the unpolarized case. Accurate distributions are now known 
from very small x to relatively large x. The same parametrization study can be done for 
the longitudinally polarized parton distributions because many experimental data are now 
available for the structure function gi [|J. 

The proton spin is so far investigated mainly in the longitudinally polarized structure 
function g\. However, it is also interesting to test it by transversely polarized structure 
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functions. In particular, the leading-twist structure function hi should be an interesting 
topic. It is expected to be measured in the transversely polarized Drell-Yan processes at 
RHIC. We should try to understand the properties of h\ before the experimental data are 
taken. The transversity distribution is an unexplored subject experimentally. If future 
experimental data on h\ are much different from theoretical expectation, a new field of spin 
physics could be developed. 

We have been working on calculation of the next-to-leading-order (NLO) anomalous 
dimensions for h\ and also on numerical solution of its evolution equation ||. We also 
studied the numerical solution of the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
evolution equations for the longitudinally polarized distributions 0. In general, parton 
distributions are provided at Q 2 ~ 0(1 GeV 2 ) and experimental data are taken at different 
Q 2 points. Therefore, our studies on the numerical solution should be useful for extracting 
information on the polarized parton distributions from polarized RHIC data. 

The major audience are traditional nuclear physicists at this symposium, so that we first 
discuss how each unpolarized parton distribution is obtained from various experimental data 
in section [| The section |3| is devoted to the activities of the RHIC-Spin-J working group. A 
brief outline of the parametrization project is explained. We discuss the numerical solution 
of the DGLAP Q 2 evolution equations for the longitudinally polarized and transversity 
distributions in section |j. The summary is given in section [5[ 

2 Extraction of parton distributions 

Nucleon substructure has been investigated through various high-energy experiments. For 
example, electron or muon deep inelastic scattering has been used. Its unpolarized cross 
section is related to two structure functions Fi and F 2 . They depend in general on two 
kinematical variables Q 2 = —q 2 and x = Q 2 /2p ■ q where q is the virtual photon momentum 
and p is the nucleon momentum. It is known that the structure functions are almost inde- 
pendent of Q 2 , which is referred to as Bjorken scaling. It indicates that the photon scatters 
on structureless objects, which are called partons. They are now identified with quarks and 
gluons. Because the variable x is the light-cone momentum fraction carried by the struck 
quark, it is possible to extract quark-momentum distributions from the F 2 data. 

The x dependence of the parton distributions are expressed by several parameters at 
fixed Q 2 {=Qq)~- for example x fi(x, Qq) = Air]iX ai (1 — (1 + ji x + pi \/x), and the 
parameters are optimized so as to explain the existing data. There are abundant experimen- 
tal data on unpolarized parton distributions such as those taken by electron/muon/neutrino 
deep inelastic scattering, Drell-Yan, and direct-photon process. With these experimental 
data, accurate parton distributions are now known although there are still detailed issues 
on the distributions at large x and very small x. The updated information on the opti- 
mum parton distributions is given at |http:/ /durpdg.dur. ac.uk/HEPDATA/HEPDATA. html| 
and |http:/ /www.phys.psu.edu"7T cteq. The major parametrization groups are CTEQ (Coordi- 
nated Theoretical/Experimental Project on QCD Phenomenology and Tests of the Standard 
Model), GRV (Gliick, Reya, and Vogt), and MRS (Martin, Roberts, and Stirling). In this 
section, we introduce how the parton distributions are extracted from various experimental 
data for a novice. 
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2.1 Valence-quark distributions 



The parton distributions are usually separated into valence-quark, sea-quark, and gluon 
distributions. The valence-quark distribution is defined by the difference between quark and 
antiquark distributions: q v = q — q. It can be obtained by analyzing neutrino deep inelastic 
data. The cross section of the neutrino reaction v^+p — > +X is expressed as da oc £^ U W^ U 
with the lepton tensor £ pu and hadron tensor W^ v J7J. The lepton tensor is calculated as 

&» = u{k)-f(l - l5 )u{k') ■ u(k') 7 "(l - 7 5 )u(k) 

spins 

= 2 [k»k' v + k u k' tl -g^k- k' + ie^ pa k p k' a ] . (2.1) 



Because of the axial vector term 7^75 in the neutrino reaction, there appears the antisym- 
metric term e pvpo ", which does not exist in the electron scattering. Due to the existence of 
this term in the lepton tensor, the hadron tensor could also have an antisymmetric one: 

ww = -Wi -*jr) + w *jp ^ - ^V) V - ^rO - ii w * w p p q a ■ (2.2) 



From these expressions, the cross section becomes 



da ± Gl E' 2 



dE'dtt 2tt 2 (1 + Q 2 /M, 



2 ^2 
w) 



2W 1 (v, q 2 ) sin 2 ^ + W 2 (v, q 2 ) cos 2 ^ =F ^^W 3 sin 2 ^ 



(2.3) 



where ± indicates W ± exchange processes. The new term with W3 is associated with the 
parity violation, and it is given by the difference between the left and right cross sections 
for the W: W 3 oc a L — a R . Instead of W3, the scaling function F 3 = z/W 3 is usually used. 
Because of the a L — a R property, the structure function F 3 is given as the quark — antiquark 
distribution form: F 3 P = 2(d + s — u — c) and F 3 P = 2{u + c — d — s). Hence, we have 

F u P , pup 

F/= 3 ^ 3 =u v + d v , (2.4) 

if s = s and c — c. It is known theoretically, for example, that the s distribution is not 
exactly equal to the s. However, because the corrections are expected to be small, the valence 
quark distribution u v + d v can be extracted from the neutrino data. 



2.2 Sea-quark distributions 

The F 2 structure function can be used for extracting both valence and sea quark distributions. 
It is given by quark distributions weighted by the square of their charges: 

F 2 (x, Q 2 ) = e i * I <l(x, Q 2 ) + fa Q 2 ) ) > (2-5) 

i 

in the leading order. Because it is dominated by the valence distributions at medium and 
large x regions, q v can be obtained directly from the F 2 data in these regions. The F 2 
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structure function at small x is , on the other hand, dominated by the sea-quark distributions. 
Therefore, q sea (= q) can be obtained in the small x region from the F 2 data. 

The other method is to use Drell-Yan processes p\ + p 2 — > + X. Its cross section 

is described by the qq annihilation processes: 

where Q 2 is the dimuon mass squared: Q 2 = fn 2 ^, and r is given by r = m 2 /s = X\x 2 . 
The rapidity y is defined by dimuon longitudinal momentum P£ and dimuon energy E* in 
the cm. system: y = (l/2)ln[(E* + P£)/(E* — PI)]. The momentum fractions X\ and £2 
can be written by these kinematical variables: X\ = \fre v and x 2 = \fre~ y '. In the large 
x F = X1—X2 region, namely at large x±, the cross section is dominated by the projectile-quark 
annihilation process with an target antiquark: da oc q v (xi)q(x2) because of q(xi) « 1. This 
equation indicates that the sea-quark distribution q(x 2 ) can be measured if the valence-quark 
distributions in the projectile are known. 

In these days, the details are investigated for the sea-quark distributions. Neutrino- 
induced dimuon events are interpreted by the charm-production process: + (s, d) — >yT + c 
then c — > s + fi + + so that the strange quark distribution can be extracted from the data. 
It is now known s + s ~ 0.5 (u + d) at Q 2 ~ 20 GeV 2 . Furthermore, flavor asymmetry in 
light antiquark distributions was also found by the Gottfried-sum-rule violation and by the 
p-n asymmetry in the Drell-Yan experiments ||. 



2.3 Gluon distribution 

There are two major methods for finding the gluon distribution. The first one is to use 
scaling violation data of the F 2 structure function, and the second is to use direct photon 
data. We mentioned that the structure functions are almost independent of Q 2 . However, it 
is known that they vary as a function of InQ 2 even though it is small. This is called scaling 
violation. The Q 2 dependence is described by the DGLAP evolution equations, which are 
the coupled integrodifferential equations: 

9 ( q s (x,Q 2 ) \ _ a s (Q 2 ) ( P qq (x,Q 2 ) P q9 (x,Q 2 ) \ ( q s (x,Q 2 ) \ . ? , 
dlnQ 2 \ g(x,Q 2 ) J ~ 2vr V PU^Q 2 ) P 99 (^Q 2 ) J V 9(x,Q 2 ) J ' [ZJ) 

where q s = J^iiVi + Qi) is the singlet quark distribution. The convolution integral ® is defined 
by f( x ) ®g{ x ) — J x f( x /y) g{y) dy/y. Each term describes the process that a parton pj with 
the nucleon's momentum fraction y splits into a parton pi with the momentum fraction x 
and another parton. The splitting function P PiPj (z) determines the probability that such 
a splitting process occurs and the p^-parton momentum is reduced by the fraction z. The 
Q 2 dependence of the quark distribution is affected by the gluon distribution through the 
splitting function. In the small x region, the right-hand sides of the evolution equations 
are dominated by the gluon terms. Therefore, the gluon distribution can be extracted from 
the Q 2 dependent data of F 2 . To be explicit, the Q 2 dependence of F 2 is approximately 
proportional to the gluon distribution: dF 2 (x, Q 2 )/d(ln Q 2 ) w 10 a s G(2x, Q 2 )/(27tt). 
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The second method is to use direct photon 
data. The leading-order cross section is calcu- 
lated by the QCD Compton scattering (gq — > 
jq) and annihilation (qq — > jg) processes. Be- 
cause the gluon distribution is involved in the 
Compton process in Fig. [I], it is possible to ex- 
tract it from the experimental data. The quark 

and antiquark distributions are also involved in 

,-i ,1 i i 11 .Figure 1: QCD Compton process, 

the reaction; however, they are relatively well 

known from deep inelastic scattering data. It is also possible to isolate the annihilation 
contribution by using both pp and pp reaction data. From the scaling violation and direct 
photon data, the gluon distribution in the nucleon is now well known except for the small 
and large x regions. The large x distribution is particularly important for judging whether 
the CDF jet production data could suggest "subquark" structure. It is desirable to measure 
the gluon distribution in the x region, x ~ 0.5, directly by some kind of experiments. On 
the other hand, the small x distribution is not well fixed as it is obvious from the Rl, R2, 
R3, and R4 type distributions in the MRS parametrizations ||. It is partly because a wide 
range of Q 2 data is not available in the structure function F2. These x regions of the gluon 
distribution should be studied experimentally in future. 



3 Activities of the RHIC-Spin-J working group 

With the completion of the RHIC facility in mind, we created a working group within Japan 
for studying possible spin physics at RHIC. We have a monthly meeting at Riken among 
related theorists and experimentalists. One of the activities is to study the parametrization 
of the polarized parton distributions. Even though the study is still preliminary, we discuss 
an outline of our efforts. 

The parametrization of the unpolarized distributions has been studied extensively as 
discussed in the previous section. On the other hand, the polarized distributions are not 
well known yet due to lack of a variety of experimental data. However, several sets of 
longitudinally polarized distributions are proposed by analyzing the g\ structure functions 
for the proton, deuteron, and 3 He |4j]. With stimulation of these studies and recent polarized 
measurements, we initiated a parametrization project within Japan [Q. Because this is a 
part of the Japanese RHIC-Spin activities, the ultimate purpose is to use our studies for the 
RHIC-Spin experiments. At this stage, the group consists of Y. Goto, N. Hayashi, M. Hirai, 
H. Horikawa, S. Kumano, M. Miyama, T. Morii, N. Saito, T.-A. Shibata, E. Taniguchi, and 
T. Yamanishi. For the time being, this collaboration intends to obtain optimum polarized 
parton distributions for explaining all the available experimental data. We have been working 
for a year, and we should be able to complete our first work in a few months. It includes 
complete NLO analyses of the gi structure functions. The obtained distributions will be used 
for experimental studies at RHIC. Furthermore, once new experimental results are obtained 
at RHIC or at other facilities, we try to reanalyze the data for getting updated distributions. 

We work in three subgroups: data analysis, parametrization, and Q 2 evolution. The data 
analysis subgroup (Y. Goto, N. Hayashi, N. Saito, T.-A. Shibata, and E. Taniguchi) collects 
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all the available experimental data and investigate their systematic uncertainties. At this 
stage, the data come from measurements of the g\ structure functions for the proton and 
"neutron". The parametrization group (T. Morii, H. Horikawa, and T. Yamanishi) tries 
to understand the meaning of obtained parameter values. They also spend much time for 
running the fitting program. The third subgroup (M. Hirai, S. Kumano, and M. Miyama) 
contributes to the Q 2 evolution program. The evolution equations are coupled integrodif- 
ferential equations with complicated splitting functions. This subgroup develops an efficient 
program for numerical solution of the evolution equations. Our parametrization studies were 
reported by N. Saito at the spring JPS meeting in 1997 ||. Recent studies were also reported 
by M. Hirai and T. Yamanishi at this meeting j3|. 

The distributions are set up at certain Q 2 (= Qq). Experimental data are taken, in 
general, at different Q 2 points. The DGLAP evolution equations are used for calculating 
the distribution variations from Qq to experimental Q 2 . The DGLAP equations assume 
perturbative QCD, so that both Qq and Q 2 have to be in the perturbative region. The 
Qq is taken typically in the 1—4 GeV 2 region except for the GRV distributions, where 
small Qq (~0.3 GeV 2 ) is chosen. If it is possible, it is desirable to take very large Qq where 
perturbative QCD is certainly valid. However, because many data are taken in the Q 2 region 
1 GeV 2 < Q 2 <a few dozen GeV 2 , we end up taking relatively small Qq at this stage. 

We provide polarized parton distributions with a number of parameters at Qq, for exam- 
ple at 1 GeV 2 . These parameters are determined by fitting the existing experimental data 
on Ai, which is given by 

a c 9i(x,Q 2 ) _ 2 2x(l + R) 

Al ~ Fl (x } Q 2 )- 9l ^ Q) F 2 (x,Q 2 ) ' (3 ' 1} 

where the Fi structure function is related to the F 2 structure function and the R function by 
R = (i*2 — 2xF\)/ (2xFi). The g\ structure function is calculated by the convolution integral 
of the polarized parton distributions with the coefficient functions: 

9i{x, Q 2 ) = \^2^ q \ f — &C q (x/y, Q 2 ) [Aq(y, Q 2 ) + Aq(y, Q 2 ) ] 

+ £ d j- AC g (x/y, Q 2 ) Ag(y, Q 2 ) j . (3.2) 

It is not a good idea to analyze the published g± values directly because it is often assumed 
in the experimental analyses that A\ is Q 2 independent. As it is obvious from the results in 
the next section (see Figs. ^ and |5|), it is not a correct assumption. However, considering 
the present experimental errors, one cannot help using the assumption. In our analysis, 
the asymmetry A\ data are analyzed without using the Q 2 independent assumption. We 
investigated various functional forms of the initial distributions; however, the studies have not 
been completed yet. For example, the x dependence is given as x Afi(x, Qq) = Ai rji x ai (1 — 
x)^ (l+7i x+pi yfx) in the analysis of Gehrmann and Stirling (GS). The F 2 structure function 
is given in the similar way as 
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F 2 (x, Q 2 ) = J2 4 x { / - C q (x/y, Q 2 ) [ q(y, Q 2 ) + q(y, Q 2 ) ] 

q I J x y 

+ fjC g (x/y,Q 2 )g(y,Q 2 )]. (3.3) 

Because the unpolarized distributions are not 
studied by our group at least at this stage , we 
take distributions fi(x, Qq) by CTEQ, GRV, or 
MRS. With these preparations, we are ready for 
the optimization. We compare the theoretical 
asymmetries A\ with the experimental data and 
minimize 



x 2 = £ , ~ i ' ( 3 - 4 ) 

by using the CERN subroutine MINUIT. Because the only information comes from the lon- 
gitudinally polarized structure function gx, the polarized sea-quark and gluon distributions 
are not well determined at this stage. We tried to find the optimum parameter set by fitting 
the data with the Q 2 evolution effects. In Fig. ^], our preliminary fitting results are shown 
with various experimental data for gf ||. The dotted curve is the initial distribution at 
Q 2 =l GeV 2 , and the solid one is the distribution at Q 2 =10 GeV 2 . We are still working; on 
obtaining the best fit. 

We have been studying the parametrization for almost a year. The leading order (LO) 
fitting program was completed a half year ago, and we try to understand physics meaning 
of the obtained parameters. On the other hand, we completed the NLO program and the 
NLO analysis is in progress. We should be able to report our results within a few months. 




4 Numerical solution of longitudinal and transversity 
Q 2 evolution equations 

The DGLAP evolution equations are frequently used in theoretical and experimental analy- 
ses. In fact, our evolution program is used in the parametrization studies in the last section. 
The evolution equations are coupled integrodifferential equations with complex splitting 
functions in the NLO case. It is, therefore, important to create a computer program to solve 
them accurately and to publish it so that other people could use the subroutine. We have 
been working on this topic for several years j5], [| 0. In this section, we report the numer- 
ical solution of the Q 2 evolution equations for the longitudinally polarized and transversity 
distributions. 

The longitudinal DGLAP equations are given in the similar way with the unpolarized 
ones. The nonsinglet equation is 
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° Aq NS (x,Q 2 ) =^AP 9± »®A^(,,Q 2 ), (4.1) 



dlnQ 2 ^ sv y 2tt 

where Aq NS = xAq NS is a longitudinally polarized nonsinglet parton distribution with Aq = 
g-f — gj,. The parton distribution and the splitting function multiplied by x are denoted 
as f(x) = xf(x). The function AP q ± NS = xAP q ± NS is the polarized nonsinglet splitting 
function. The notation q^ 1 in the splitting function indicates a Aq + = Aq+Aq or Aq~ = Aq— 
Aq distribution type. The singlet Q 2 evolution is described by the coupled integrodifferential 
equations due to the gluon participation: 



d ( Aq s (x,Q 2 ) \ _ a s {Q 2 ) ( AP qq (x,Q 2 ) AP qg (x,Q 2 ) \ / Aq s (x,Q 



)2\ 



® A~?n2 -(4-2 



«91nQ 2 V A^?(^Q 2 ) J ~ 2tt \ AP gq (x, Q 2 ) AP gg (x,Q 2 ) J \ Ag(x,Q 

The singlet quark distribution is defined by Aq s = x Aqf where i is the flavor, and 
Ag(x, Q 2 ) is the gluon distribution A'g = x (g + \ — g-i). The LO and NLO Jll| evolution 
equations are described by the above integrodifferential equations. The differences between 
the LO and NLO are contained in the running coupling constant a s and in the splitting 
functions. The renormalization scheme is the modified minimal subtraction scheme (MS) 
in the NLO analysis. 

In addition to the unpolarized and longitudinally polarized evolution equations, the NLO 
evolution for the transversity distribution was completed recently ||12|| . The transversity 



distribution is denoted as hi, Sq, or A T q. Throughout this paper, the notation A T q is used. 
Because the gluon does not couple due to the chiral-odd nature of A T q, the evolution is 
described by a single integrodifferential equation, 

(> A T g ± (x,Q 2 ) ='^l^ T P q± (x)®A T q ± {x,Q 2 ) . (4.3) 



dlnQ 2 " ^ ' 2tt 

Numerical analysis of the transversity evolution is simple in the sense that it is given only 
by the "nonsinglet" type equation. 

The longitudinally polarized and transversity distributions will be measured at RHIC by 
polarized p+p reactions. However, they are measured at different Q 2 points. In order to 
extract information on the polarized parton distributions, an accurate program for solving 
the evolution equations should be developed. We have been studying on this topic, and the 
following discussions are based on the numerical results in Refs. 0, |6| . We simply divide the 
variables Q 2 and x into small steps. Then, the integration and differentiation are defined by 

df(x) = f(x m +i) - f(x m ) ^ 4 4 ^ 

N x 

f(x) d,X=^2 5X ™ f( X m) ■ (4.5) 
m=l 

With these replacements, the evolution equations could be solved rather easily. We call 
this method a "brute- force" method M. The better method, by means of convergence and 
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Figure 5: Spin asymmetry at x = 0.25 H]. 

D 

In Fig. [|, our evolution curves at x=0.035 are shown with the asymmetry A\ data. The 
dashed and solid curves indicate the LO and NLO evolution results, respectively. In the 
large Q 2 region, both results are almost the same; however, they differ significantly at small 
Q 2 particularly in the region Q 2 < 2 GeV 2 . In the medium x region (x=0.25), the difference 
becomes larger again at small Q 2 as shown in Fig. |5[ From these figures, we find that the 
asymmetry has Q 2 dependence although it is not large. People used to assume that the 
asymmetry is independent of Q 2 by neglecting the Q 2 evolution difference between g% and 
Fx in analyzing the experimental data. We find clearly that it is not the case. For a precise 
analysis, the Q 2 dependence in the asymmetry has to be taken into account. 

Next, we discuss the transversity evolution A T q. It has a chiral-odd property so that it 
cannot be found in inclusive deep inelastic electron scattering. It is expected to be measured 
in the RHIC-Spin project by the transversely polarized Drell-Yan processes. The distribution 
A T q is the probability to find a quark with spin polarized along the transverse spin of a 
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polarized proton minus the probability to find it polarized oppositely. There is a problem in 
studying the transversity evolution in the sense that the input distribution is not available. 
However, it is known within quark models that the transversity distributions are almost the 
same as the corresponding longitudinally polarized distributions. Therefore, we may use a 
longitudinal distribution as an input transversity distribution at small Q 2 , where the quark 
models would be valid. This consideration is fine at Q 2 ~ ^qcd- O n the other hand, the 
perturbative QCD cannot be used at such a small Q 2 . In the following numerical analysis, we 
simply assume Aq = A T q at relatively large Q 2 (=4 GeV 2 ) by taking the GS-A distributions 
as the input. 

We show the evolution results of two distributions, the singlet distribution x^2 i (A T q + 
A T q) and the valence distribution x(A T u v + A T d v ). First, the singlet evolution is shown 
in Fig. []. The same GS-A distribution is assumed for the transversity and longitudinally- 
polarized parton distributions at Q 2 =4 GeV 2 . The initial distribution is shown by the 
dotted curve. It is evolved to the distributions at Q 2 =200 GeV 2 by the transverse or longi- 
tudinal evolution equation. The transversity NLO effects increase the evolved distribution 
at medium-large x and also at small x (<0.01), and they decrease the distribution in the 
intermediate x region (0.01 < x < 0.1). The transversity evolution differs significantly from 
the longitudinal one either in the LO case or in the NLO. The evolved transversity distribu- 
tion A T q s is significantly smaller than the longitudinal one Aq s in the region x ~ 0.1. The 
magnitude of A T q s itself is also smaller than that of Aq s at very small x (< 0.07). 

Second, the evolution of the valence-quark distribution is shown in Fig. [7]. The evolution 
calculations are done in the same way. The dotted curve is the initial distribution, and the 
evolved distributions are solid and dashed curves. The evolved transversity distributions are 
significantly smaller than the longitudinally polarized ones, in particular at small x. 




Figure 6: The singlet transversity evolution re- 
sults are compared with those of longitudinally 
polarized singlet distribution |5| . 



Figure 7: Evolution results of the valence- 
quark transversity distribution x(A T u v + A T d v ) 
are compared with those of the longitudinally po- 
larized one |5| . 



We have created the useful computer programs, BFP1 and H1EVOL, for solving the 
DGLAP evolution equations for the longitudinally polarized parton distributions and for 
the transversity distributions. They should be convenient for theoretical and experimental 
researchers, who are involved in high-energy spin physics. Our program can be obtained from 
M. Hirai, SK, or M. Miyama upon email request. For the details, the reader may look at the 
WWW home page: [http: / / www.cc.saga-u.ac.jp/saga-u/riko/physics/ quantuml / program.html 
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5 Summary 



Since the Riken decided to join the RHIC-Spin project, we have been studying the physics as- 
pects of the project. In particular, the parametrization of polarized parton distributions has 
been investigated among the members. For general audience at this symposium, the method 
of extracting unpolarized parton distributions is first discussed. Then, our parametrization 
project was explained although the results are still preliminary. The optimum parton dis- 
tributions are obtained by fitting existing experimental data of the spin asymmetry Ax. On 
the other hand, we have been working on creating efficient computer programs of solving 
the Q 2 evolution equations. Now, the evolution programs for the unpolarized, longitudinally 
polarized, and transversity parton distributions are available. Our efforts should be useful 
for researchers in the area of high-energy spin physics. 
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